Neutrino mixing as a source of dark energy 
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We show that the vacuum condensate due to neutrino mixing in quantum field theory (QFT) 
contributes to the dark energy budget of the universe which gives rise to the accelerated behavior 
of cosmic flow. The mechanism of neutrino mixing is therefore a possible candidate to contribute 
to the cosmological dark energy. 
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Data coming from cosmic microwave background ra- 
diation (CMBR) [1, 2], large scale structure [1, Q and 
type la supernovae [5|, used as standard candles, inde- 
pendently support the picture that the today observed 
universe can be consistently described as an accelerating 
Hubble fluid where the contribution of dark energy com- 
ponent to the total matter-energy density is CIa — 0.7. 
The big challenge is then the one of explaining such a 
bulk of dark energy component. 

On the other hand, in recent years great attention has 
been devoted to the neutrino mixing phenomenon. The 
theoretical investigation of the neutrino mixing, firstly 
propose d by Pontecorvo has been pursed in depth 
0) B S [13, Ell ) and more recently the issue of the 
construction of the flavored space of states has been set- 
tled in the framework of t he q uantum field theory (QFT) 
formalism [H, El [H, [H, El El, [11 with the discovery 
of the unitary inequivalence between the flavored vacuum 
and the massive neutrino vacuum (isl . E3 1 the associated 
finding of the neutrino-antineutrino pair condensate con- 
tributing to the vacuum energy and the new oscilla- 
tion formulas [3, El) E^ E3 ■ The recent experimental 
achievements proving neutrino oscillations f 2ll . [2^ and 
the progresses in the QFT theoretical understanding [ll 
of the neutrino mixing thus provide a challenging and 
promising path beyond the Standard Model of electro- 
weak interaction for elementary particles. 

In this paper we show that these two most interesting 
issues are intimately bound together in such a way that 
one of them, namely the neutrino mixing phenomenon, 
appears to provide a contribution, till now unsuspected, 
to the vacuum dark energy component. The structure of 
the fiavor vacuum and its unitary inequivalence to the 
vacuum for the massive neutrinos play a significant role 
in obtaining the result we report below. 

We want to stress that in this paper we do not solve 
the problem of the arbitrariness of the momentum cut- 
off. The cut-off problem should therefore not be confused 
with the main message we want to present in this paper, 
namely that the very same phenomenon of the neutrino 
mixing appears to provide a contribution to the cosmo- 
logical dark energy component. As shown below, such 
a feature, unsuspected till now, clearly manifests itself 
provided the mathematically correct QFT treatment of 



neutrino mixing [H, El El, El E3, El, El, HI is consid- 
ered. 

In the simplest explanation, the so called KCDM 
model, the cosmological constant contributes for almost 
70% to the total matter-energy density budget. The 
standard theory of cosmological constant is based on the 
fact that the vacuum zero point energy cannot violate 
the Lorentz invariance of the vacuum and therefore the 
corresponding energy-momentum tensor density has the 

form = (o|r^.|o) 

— PA5/ii/, where pA is a constant, 
i.e. a Lorentz scalar quantity. In the traditional pic- 
ture the vacuum itself can be thought as a perfect fluid, 
source of the Einstein field equations and one derives [l^ 
the equation of state pA — wp\, with p\ denoting the 
vacuum pressure and the adiabatic index w equals — 1. 
As well known [24] , one of the central pillars of Lorentz 
invariant local QFT is the very same definition of the 
vacuum state according to which it is the zero eigen- 
value eigenstate of the normal ordered energy, momen- 
tum and angular momentum operators. Therefore, ex- 
cluding by normal ordering zero-point contributions, any 
non-vanishing vacuum expectation value of one of these 
operators signals the breakdown of Lorentz invariance, 
since the vacuum would be dependent on space and/or 
time. The Lorentz invariance vacuum therefore implies 
^uac = : : |0) = 0, (as usual normal ordering is 
denoted by the colon : ... :). 

Usually, in the jargon one roughly expresses the 
Lorentz invariant characterization of the QFT vacuum, 
by saying that preserving the Lorentz invariance re- 
quires to exclude that kinematical terms in the energy- 
momentum tensor may contribute to the vacuum expec- 
tation values. 

Suppose that, as said above, the contribution of the 
(zero point) vacuum energy density is taken to be equiv- 
alent to that of the cosmological constant A, which is 
expressed by pA — A/{8ttG). Then, however, it turns 
out that the vacuum expectation value of the energy- 
momentum tensor is divergent, both for bosonic and 
fcrmionic fields, and this shortcoming can be addressed as 
the cosmological constant problem. By choosing to regu- 
larize the energy-momentum tensor by an ultraviolet cut- 
off at Planck scale, one gets a huge value for the vacuum 
energy density pyac ^ c^/G^h - lO^^Gey" which is 123 
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orders of magnitude larger than the currently observed 
Pa ~ 10~^'^GeV'^ [25|]. Also using a quantum chromo- 
dynamics (QCD) cut-off [1^ the problem is not solved 
since p^*^^ ~ 10~^GeV'^ is still enormous with respect 
to the actual observed value. 

Furthermore, there is another aspect which has to be 
taken into account: observations point out that cosmic 
flow is "today" accelerating while it was not so at inter- 
mediate redshift z (e.g. 1 < z < 10). This situation 
gave rise to structure formation during the matter dom- 
inated era [271, 28]. This is an indication of the fact that 
any realistic cosmological model should roughly undergo 
four phases: an early accelerated phase (inflation), inter- 
mediate decelerated phases (radiation and matter dom- 
inated) and a final, today observed, accelerated phase. 
Obviously, the dynamical evolution (time dependence) of 
the cosmological constant through the different phases, 
breaks the Lorentz invariance of the vacuum and one has 
to face the problem of the dynamics for the vacuum en- 
ergy in order to match the observations. In this case, 
therefore, we are not properly dealing with cosmological 
constant. Rather, we have to take into account some form 
of dark energy which evolves from early epochs inducing 
the today observed acceleration. Such a dynamical evo- 
lution of the dark energy, namely time dependence of the 
energy vacuum expectation value, violates the Lorentz 
invariance. 

In the literature there are many proposals to achieve 
cosmological models justifying such a dark energy compo- 
nent, ranging from quintessence |29], to braneworld 
to extended theories of gravity [31]. These approaches 
essentially consist in adding new ingredients to the dy- 
namics (e.g. scalar fields), or in modifying cosmological 
equations (e.g. introducing higher order curvature terms 
in the effective gravitational action). 

In this letter we show that, due to the condensate 
of neutrino-antineutrino pairs, the vacuum expectation 
value of the energy-momentum tensor naturally provides 
a contribution to the dark energy p™^, which in the 
early universe satisfies the strong energy condition (SEC) 
Pvac + ^Pvac — 0, and at present epoch behaves approx- 
imatively as a cosmological constant. Here is the 

vacuum pressure induced by the neutrino mixing. Un- 
der such a new perspective, the energy content of the 
vacuum condensate could be substantially interpreted as 
dynamically evolving dark energy. 

The main features of the QFT formalism for the neu- 
trino mixing are summarized as follows. For the sake 
of simplicity, we restrict ourselves to the two flavor case 
[3] ■ Extension to three flavors can be found in Ref. [13] ■ 
The relation between the Dirac flavored neutrino fields 

(x) , i^fj, {x) and the Dirac massive neutrino fields vi (x) , 
h'2{x) is given by 

f Veix) \ _ ( cosd sme\ ( vi{x) \ , . 

\ v^,{x) ) " \-sine cose) \ V2{x) ) 

being 9 the mixing angle. The mixing transformation 
([T|) can be written as V(r{x) = Gg^{t) Vi{x) Gg{t), 



where (cr, i) = (e, 1), (/x, 2), and Ge{t) is the transfor- 
mation generator. The flavor annihilation operators are 
defined as al^^{t) = Gg\t) a^^ Gg{t) and /31k,^(t) = 
Ge'^(0 /?!lk,j Geit). They annihilate the flavor vacuum 
\Q{t))e,ii = Gg^{t) \0)i,2, where l0)i^2 is the vacuum an- 
nihilated by j and /3^j^ 

The crucial point of our discussion is that ]0(t))e,^, 
which is the physical vacuum where neutrino oscillations 
are experimentally observed, is a (coherent) conden- 
sate of Qfk.i (/3k, i) neutrinos (antincutrinos): 

(2) 

where i = 1,2, the reference frame k = (0,0, ]k]) has 
been adopted for convenience, is the Bogoliubov coef- 
ficient enterin g th e mixing transformation (see for exam- 
ple Refs. P^TitI) and \0)e4i denotes \0{t))e^^ at a con- 
ventionally chosen time t — 0. As a consequence of its 
condensate structure the physical vacuum \0{t))e,fj. turns 
out to be unitary inequivalent to l0)i_2 [13] ■ For brevity, 
we omit here to reproduce the explicit expression of Vk 
which can be found, e.g., in Refs. [H, [i3,[ii. We only 
recall that Vk is zero for mi = m2, it has a maximum 
at ]k] = ^m^m^ and, for ]k] 3> , it goes like 

jVkl^ — ("T-2 ~ "^1)^/(4]^]^). For massless neutrinos, as 
well known, one does not have mixing (mi = m2 = 0). 
The oscillation formulas for the fiavor charges Qe,fj.{t) 
are obtained by computin g th eir expectation values in 
the physical vacuum jO)^ [ItI. ITqI]. 

Let us now calculate the contribution p™" of the neu- 
trino mixing to the vacuum energy density. We con- 
sider the Minkowski metric (therefore we use the nota- 
tion 77'*" instead of g'^'^)- The particle mixing and oscil- 
lations in curved background will be analyzed in a sep- 
arate paper. Eq. ([2]) suggests that the energy content 
of the physical vacuum gets contributions from the ak,i 
and /?_k.i neutrino condensate. Therefore, as custom- 
ary in such circumstances, we compute the total energy 
Too = / (PxTqo{x) for the fields z/i and 1/2, 

■ ■■-T.f ^'kc..,. (<,aL, + , (3) 

r '' 

with z = 1, 2 and where : ... : denotes the normal ordering 
of the Qfk.i and /3_k.i operators. Note that T™ is time 
independent. 

We remark that we have e^^(Ol : T^™ : lO)^,^ = 
e..t.{0{t)\ ■■ : |0(t))e,M, for any t, within the QFT 

formalism for neutrino mixing. The contribution p,™" 
of the neutrino mixing to the vacuum energy density is 
thus obtained: 

ie,^(Ol^:T™:lO)e,p=p;:?r?o°. (4) 

i 

By using Eq.®, we then have 
2 

p™^ = -sin^^ / dkk\cok,i+u;k,2)\V^\^ (5) 

TT Jo 
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where the choice of the cut-off K will be discussed below. 
Similarly, the expectation value of T^^i^ in the vacuum 

|0)e.^ gives the contribution p^^^ of the neutrino mixing 
to the vacuum pressure: 



IE 



rpJJ . |n\ _ ^mix 



(6) 



where no summation on the index j is intended. Being, 
for each diagonal component, 



(no summation on repeated indices), in the case of the 
isotropy of the momenta: = k'^ — k^, 3(fc-')^ — fc^, we 
have T^^ = T^^ = T^'^ and the following equation holds 



Pv 



37r 



K 



dkk" 
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l^k|'.(8) 



Eqs.® and ([5]) show that Lorentz invariance is broken 
and jO™^ ^ —Pmc foi' value of the masses mi and 
m2 and independently of the choice of the cut-off. We 
observe that w ~ —1/3 when the cut-off is chosen to 
he K TOi,m2, cf. Eqs.® and ([8]) and the discussion 
below for the choice of K. 

It is worth stressing that the violation of the Lorentz 
invariance originates from the neutrino-antineutrino con- 
densate structure of the vacuum. Indeed, as it appears 
from the computation reported above, in the absence of 
such a condensate, i.e. with |Vkp = 0, the vacuum ex- 
pectation value of each of the (0, 0) and (j, j ) components 
of the energy-momentum tensor would be zero. We also 
remark that the non-zero expectation value we obtain is 
time-independent since, for simplicity, we are consider- 
ing the Minkowski metric. When the curved background 
metric is considered, |VkP gets a dependence on time, 
as we will show in a forthcoming paper. In any case, 
the contribution to the vacuum expectation value of T'^'^ 
is found to be non- vanishing, in the present computa- 
tion (or in the curved background case), not because the 
adopted metric is flat (or not), but because of the non- 
trivial structure of physical vacuum due to the mixing 
phenomenon (which manifests itself in the non- vanishing 
of |FkP). 

The above result holds in the early universe, when the 
universe curvature radius is comparable with the oscilla- 
tion length. At the present epoch, in which the breaking 
of the Lorentz invariance is negligible, the non-vanishing 
vacuum energy density p™" compatible with Lorentz in- 
variance cannot come from condensate contributions car- 



rying a non-vanishing 9^ 



kn 



{ujk ,kj), as it happens 



in Eqs.© and © (see also Eqs.® and (O). This means 
that it can only be imputed to the lowest energy contri- 
bution of the vacuum condensate, approximatively equal 
to 



mix 

Pa 



3 e,M\ 



Vi{x)vi{x) : |0)e,^. (9) 



Consistently with Lorentz invariance, the state equation 
is now p™" ~ —Pa", where explicitly 
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The result pI7|) shows that, at present epoch, the vac- 
uum condensate, coming from the neutrino mixing, can 
contribute to the dark energy component of the universe, 
with a behavior similar to that of the cosmological con- 
stant [13 . 

We observe that, since, at present epoch, the charac- 
teristic oscillation length of the neutrino is much smaller 
than the radius of curvature of the universe, the mixing 
treatment in the flat space-time, in such an epoch, is a 
good approximation of that in FRW space-time. More 
interesting is also the fact that, at present epoch, the 
space-time dependent condensate contributions, carry- 
ing a non- vanishing fc^, are missing (they do not con- 
tribute to the energy-momentum vacuum expectation 
value). The modes associated to these missing contribu- 
tions are not long-wave-length modes and therefore they 
are negligible in the present flat universe, i.e with respect 
to the scale implied by an infinite curvature radius. 

Eqs. ([5]) and pO]) show that the contribution to the 
dark energy induced from the neutrino mixing of course 
goes to zero in the no-mixing limit, i.e. when the mixing 
angle 9 = Q and/or mi = m2, and clearly also for mass- 
less neutrinos. However, those equations also show that 
the contribution depends on the specific QFT nature of 
the mixing: indeed, it is absent in the quantum mechan- 
ical (Pontecorvo) treatment, where Vk is anyhow zero. 
This confirms that the contribution discussed above is a 
genuine QFT non-perturbative feature and it is thus of 
different origin with respect to the ordinary vacuum en- 
ergy contribution of massive spinor fields arising from a 
radiative correction at some perturbative order [32i] . This 
leads us to believe that a neutrino-antineutrino asym- 
metry, if any, related with lepton number violation [33|, 
would not affect much our result. We will consider the 
problem of such an asymmetry in a future work. 

We call the reader attention on the fact that the dis- 
cussion and the related results presented till now, which 
constitute the core of our message, are independent of 
the choice of the cut-off K . The cut-off problem is a dis- 
tinct problem to which we do not have a solution at the 
present. 

As shown in Ref . 33] , in a dense background of neu- 
trinos, as in the case of the early universe during the Big 
Bang Nucleosyntesis, flavor particle-antiparticle pairs are 
produced by mixing and oscillations with typical momen- 
tum k ~ mi+m.2 ^ ^]_^g average mass of the neutrinos. In 
this connection, we note that in the literature [IS, 113] it 
has been noticed that the observed order of magnitude 
pTnix ^ lO~*'^GeV'^ can be reproduced by cutting the 
momentum range at the QFT neutrino scales y/mim2 or 
i 0.3, nii of order of lO'^^eV, so that 



with sin- 
Sm? = mi — m? 



— mf ~ 8 X 10 °eV . This reflects the pres- 
ence of the |Vkp factor, with its momentum dependence. 
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in the above integrations, which points to the relevance 
of soft momentum (long-wave-length) modes. 

In conclusion, we have shown that the vacuum con- 
densate due to neutrino mixing contributes to the dark 
energy budget of the universe. Different behaviors of the 
vacuum expectation value of the energy-momentum ten- 
sor have been discussed referring to different boundary 



conditions in different universe epochs. Our result is in- 
dependent of the choice of the cut-off. 
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